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ABSTRACT 

We study the luminosity function of satellite galaxies around isolated primaries using 
the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy samples. 
We select isolated primaries from the spectroscopic sample and search for potential 
satellites in the much deeper photometric sample. For primaries of similar luminosity 
to the Milky Way and M31, we are able to stack as many as ~ 20, 000 galaxy systems to 
obtain robust statistical results. We derive the satellite luminosity function extending 
almost 8 magnitudes fainter than the primary galaxy. We also determine how the 
satellite luminosity function varies with the luminosity, colour and concentration of 
the primary. We find that, in the mean, isolated primaries of comparable luminosity 
to the Milky Way and M31 contain about a factor of two fewer satellites brighter than 
My = -14 than the average of the Milky Way and M31. 

Key words: Galaxies: dwarf, Galaxies: structure, Galaxies: luminosity function, mass 
function, Galaxies: Local Group, Galaxies: fundamental parameters 
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1 INTRODUCTION 

The ACDM model predicts that structure forms in a hier- 
archical manner. Large spiral galaxies like the Milky Way 
(MW) and M31 form within extended dark matter halos 
from the merging and accretion of smaller subhalos. Smaller 
structures falling into bigger haloes can survive there as sub- 
structures and host observed satellite galaxies. A strong pre- 
diction of the theory, borne out by high resolution N-body 
simulations, is that a very large number of suc h dark matter 
substructures should survive in galactic halos |Klypin et al.l 
19991: iMoore et al. hl999l : iDiemand. Kuhlen fc Madaull2007l ; 
Springel et al.ll20*0£ ). Empirical tests of this prediction have 
so far been restricted to a single system, the Local Group, 
the only one for which an estimate of the sate llite luminos- 
ity fu nction (LF) is read il y ava ilable. Indeed, iKlypin et al.l 
1 19991 ) and IMoore et all \ 19991 ) noted that the observed 
number of satellites around the MW and M31 is much 
smaller than the number of predicted substructures, giv- 
ing rise to the so-called "missing satellites problem" of the 
ACDM model. 

In the past decade or so, fainter satellites around 
the MW and M31 have been discovered in the SPSS 
(e.g. iBelokurov et alJl2008l. l20ld: |Grebelll2000l: I Irwin et al.l 



2007l:lLiu et al.ll2008l:lMartin et al]|200d:lMartin et al.ll200S 
Simon fc Gehal 120071; Ivan den Berghl l200d ; IWatkins et al] 
20091 : IZucker et all |2004 l200rj . 120071 ), but the number is 
still orders of magnitude smaller than the predicted num- 
ber of surviving cold dark matter subhalos. A number of 
solutions to this problem have been proposed. Some in- 
voke a different kind of dark matter, warm dark mat- 



ter, in which case the number of surviving substructures 
is dramatically reduced by a cut o ff in the primordial 
powe r spectrum (Moore et alj 2000l: Spergel fc Steinhardtl 



200d: lYoshida et al l l200d; iBode Ostriker, fc Turokl l200ll : 



Craig fc Davisll200ll ; lLovell et al.ll2011f ). Others retain cold 



dark matter and appeal to galaxy formation processes, such 
as photoionization and supernova feedback, to inhibit star 
formation in small subhalos thus rendering most of them 
inv isible. This idea, first m entioned nearly 20 years ago 
by iKauffmann et all l|l993h . was worked out in detail a 
decade l ater using analytical arguments and semi-analytical 
models dBullock, Kravtsov fc Weinberg! l200d : H enson et al.l 
120021 : ISomervillell2002l )~ 

The discovery of new Local Group satellites in the 
SDSS has stimula ted further obse r vatio nal an d, particularly, 
theore tical work. iKoposov et al. I (|2008T l and iTollerud et al.l 
(2008) extended the estimate of the satellite LF of the 
MW and M31 to faint magnitudes, accounting for the 
survey magnitude limit and modelling the radial den- 
sity profile of the satellite distribution. This extension 
to faint magnitudes agrees remarkably well w ith the 
ACDM model predictions of iBenson et al.l (|2002h . a re- 
sult that has been confirmed in rece nt work using related 
semi-analytic mod e lling techniques (Koposov et all 1 2009 1: 
Mufioz et alj |2009| ; iBusha et alj l20ld: ICooper et al.l hold; 
Maccio et aLlhoiol ; \Ll De Lucia, fc Helmill20ld : iFont et al l 
201lT ). Full N-body /gasdynamic simulations have also 



been car ried out to inv e stigat e the physics of satellite 
galaxies (jLibeskind et~ai1 120071: [Okamoto fc Frenk I 120091 : 
lOkamotoet al.ll20ld : IWadepuhl fc Springe]||20ld ) although 
currently these simulations only resolve the brightest exam- 
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pies. In spite of this broad agreement, interest ing discrepan- 
cies ex ist. For example, the original model of I Benson et al. 
2002), as well as the more recent model by iGuo et al.l 



2010) , rarely produce satellites a s bright as the LMC and 



SMC jBovlan-Kolchin et alJlioiol ). 

The large body of work on satellite galaxies reflects the 
importance of these objects as a critical test of the ACDM 
model on small scales. Yet, all conclusions to date regarding 
the validity or otherwise of the model rely on comparison 
with data for a few dozen satellites around just two galax- 
ies, the MW and M31. There is no guarantee that these are 
typical and indeed there is good evidence that the satel- 
lites of the two galaxies have different structural properties 
l|McConnachie fc Irwinl 120061 ; ICollins et all |2010| ). Clearly, 
robust and reliable tests of cosmological and galaxy forma- 
tion models require comparison with statistically represen- 
tative samples of galaxies and their satellites. 

Analyzing the satellite sytems of external galaxies 
is challenging because typically only o ne or two satel- 
lites are detected per primary galaxy Jrlolmberel Il969l : 
lLorrimer et~aH 11994 ; IZaritskv etabl Il993l . Il997bl ). In addi- 
tion, the real space position of the satellite with respect 
to its primary is uncertain. To circumvent the first prob- 
lem, these authors developed the method of stacking the 
primaries and their satellites in order to obtain a fair and 
complete sample which can yield statistically robust results 
for certain classes of primary galaxies. These early studies 
were limited by the relatively small samples available at the 
time. With the advent of large galaxy redshift surveys such 



as th e 2dF Galaxy Redshift Survey (2dFGRS; Collcss ct al 



2001) and the Sloan Digital Sky Survey (SDSS; IYork et al 



2000), it is now possible to construct external galaxy sam- 



ples covering a much larger volume. Studies with signifi- 
cantly improved st atistics have been carr i ed out using these 
new surveys (e.g. Sales fc Lambasl |2004| ; lYang et al.l l200d : 
lAgustsson fe Brainerdll2010h , 

In this work, we are interested in the satellite lumi- 
nosity function of specific types of isolated primary galaxies 
and, for this, the new spectroscopic surveys are still not deep 
enough. For example, even w ithin the largest galaxy redshift 
catalogue from SDSS DR7 (|Abazaiian et al.l l2009| ). where 
there are about ~ 660 000 galaxies with mj! m < 17.77, only 
a relatively small number of isolated low redshift galaxy sys- 
tems have enough de tected satellites for our purposes (e.g. 
iHwang fe Park! 20101 ) . On the other hand, the photomet- 
ric catalogue from SDSS DR7 contains ~ 96 000 000 galaxies 
with magnitudes in the m, g, r, i, z bands (roughly mj, lm < 22) 
and photometric redshifts. In this study, we used both the 
spectroscopic and photometric SDSS DR7 catalogues. To 
ensure completeness, we restrict the photometric sample in 
our main analysis to galaxies brighter than m r = 20.5 (see 
Section 4). The resulting catalogues enable us to analyze a 
sufficiently large statistical sample of galaxy systems. We 
construct o ur sample usin g meth ods similar to those de- 
veloped by lLorrimer et al.l |l994) but modified slightly to 

include photometric redshifts. 

As this project was nearing completion, iLiu et all 

<|2010l) published an investigation using similar methods 
to quantify the frequency with which satellites as bright 
as the LMC and SMC occu r around primaries similar to 
the M W. Shortly afterwards, lLares. Lambas. fe Dominguej 
i|2010l ) also published a similar study, focused on primaries 



brighter than M r < —20.5, investigating how the satellite 
LF and projected density profile depend on the primary lu- 
minosity and colour. Our work complements these studies 
by including a wider range of primary luminosities and ex- 
ploring how the satellite LF depends on the properties of 
the pri mary. Also, we adopt stricter isolation criteria than 
those of lLares. Lambas. fc Dominguej (|2010T ) . We compare 
our results with those of these studies in the discussion in 
Section [S] 

The remainder of this paper is organised as follows. In 
Section [2] we describe the selection of primary galaxies and 
their satellites; in Section [3] we develop the method of esti- 
mating the satellite LF; in Section!?] we present our estimate 
of the satellite LF for different types of primary galaxy. We 
conclude, in Section [S] with a summary and discussion of 
our results. 



2 DATA AND SAMPLE SELECTION 

We build two different catalogues for our study: a smaller 
one of galaxies with spectroscopic redshifts from which we 
select the primary galaxies (hereafter the spectroscopic cat- 
alogue) and a larger one of galaxies with photometric red- 
shifts and magnitudes from which we select the neighbouring 
galaxies (hereafter the photometric catalogue) . The spectro- 
scopic catalogue is constructed from the SDSS DR7 spec- 
troscopic subsample (north galactic cap) including all ob- 
jects with high quality redshifts (zconf > 0.7 and specClass 
= 2) and a Petrosian magnitude r ^ 17.77. The photo- 
metric galaxy catalogue is from the SDSS DR7 photometric 
subsample (north galactic cap) and includes only objects 
that have photometric redshifts, none of the flags BRIGHT, 
SATURATED, or SATUR.CENTER set and model mag- 
nitudes r ^ 22.0. We select only objects with correspond- 
ing entries in the SDSS database PhotoZ table, which nat- 
urally selects galaxies and excludes stars. As galaxies with 
r ^ 17.77 are included in both SDSS catalogues, a small frac- 
tion of the photometric catalogue galaxies also have spec- 
troscopic redshifts. We use de-reddened model ugriz magni- 
tu des and k-correct all gala xies to z — with the IDL code 
of lBlanton fc Roweisl l|2007h . In addition, we also include V- 
band magnitudes estimated from g and r-band magn itudes 
assuming V = p-0.55(g-r)-0.03 (|Smith et al.ll2002l ). This 
allows us directly to compare our results with observations 
of the MW. 

For our statistical analysis, the sample of primary galax- 
ies should be not only homogeneous but also isolated. To this 
end, we adopt a series of selection criteria summarised in the 
flow chart shown in Fig.[T] First, from the spectroscopic sam- 
ple, we select primary galaxy candidates of absolute magni- 
tude, M p , in the range M c - AM bin < M p ^ M c + AM bin . 
We then reject from this list those candidates that have, 
or could have, bright neighbours whose own satellite sys- 
tem could overlap with that of the candidate. We achieve 
this by rejecting candidates that have a neighbouring galaxy 
within a projected distance of 2_Ri nnor that is brighter than 
M p + AMf a i n t, unless that neighbouring galaxy is at a 
substantially different redshift. For neighbours with spec- 
troscopic redshifts, z s , the required redshift separation is 
|^ p — j3 s > Az s , while for those with only photometric red- 
shifts, Zphot, we require \z p — z p hot| > ce p a*. Here a* is the 
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A: Is the primary galaxy magnitude, M p , in 
the range M c - AM bin < M p < M c + AM bin ? 




yes 



Reject 
central 
galaxy 



Neighbouring galaxies are from the photo- 
metric catalogue with photometric redshifts 
and possibly also spectroscopic redshifts 

B: Are there any neighbouring galaxies 
within a projected radius of 2i? inner that are 
brighter than M p + AM faint ? 

C: For each neighbour within 2_Rj nner and 
brighter than M p + AM la i nt , if it has a spec- 
troscopic redshift, is it within Az s of the cen- 
tral galaxy? 

D: For the remaining neighbours within 
2-Rinnen brighter than M p + AMf a j nt and 
without spectroscopic redshifts, do their pho- 
tometric redshifts satisfy \z p — z c \ < a p a*7 



Figure 1. A flow chart detailing the selection criteria for isolated primary galaxies. 



photometric redshift error that we adopt (see Section[3]) and 
q p is a tolerance, which we will vary. The isolation criteria 
guarantee that there are no luminous neighbouring galaxies 
that are projected within 2_Rj Ilnor of the primary, unless these 
luminous neighbours are sufficiently far away from the pri- 
mary and appear here due to a chance projection. Using the 
photometric redshift information to identify and remove true 
background and foreground galaxies significantly increases 
the number of primary galaxies retained in our sample and 
reduces the background contamination. 

After having filtered by these criteria, the remain- 
ing isolated galaxies comprise the primary galaxy cata- 
logue. We briefly summarise the properties of this cata- 
logue. The number of primary galaxies not only depends 
on their absolute magnitude, but also on the isolation pa- 
rameters. The stricter the isolation criteria we take, the 
fewer primary galaxies we have. In the V-band, with a pa- 
rameter Set {M c , i?inncr, -Router, AA/ bin , AM faint , Az s , Q p } = 

{-21.0,0.3 Mpc,0.6 Mpc, 0.5, 0.5, 0.002,2.5]Q, the number 
of candidates is 202 351, which, after applying the isolation 
criteria, is reduced to 21 346 or about 10% of the galaxies in 



The parameter, -Router, is defined below 



Table 1. Properties of the primary galaxy samples for the 
following default choices of values for the sample selection 
parameters, {-A/ c , i?i nnC r; Router , AM bin , AM faint > A2 a , Op } = 
{M c ,0.3 Mpc, 0.6 Mpc, 0.5,0.5, 0.002, 2.5}. Quantities listed for 
each bin of V-band absolute magnitude M c are: the number of 
primary galaxy candidates (galaxies within the absolute magni- 
tude bin), the number of primary galaxies that pass all the isola- 
tion criteria, their median redshift and redshift range. 



Mc 


primary 


primaries 


median 


redshift 








candidates 




redshift 


range 






-19.0 


35893 


88 


0.043 


0.021 < 


z 


< 0.066 


-20.0 


104907 


2661 


0.105 


0.020 < 


z 


< 0.068 


-21.0 


202351 


21346 


0.098 


0.016 < 


z 


< 0.164 


-22.0 


94287 


51733 


0.142 


0.022 < 


z 


< 0.391 


-23.0 


51686 


26982 


0.203 


0.031 < 


z 


< 0.522 



this magnitude bin. The primary galaxy redshifts lie in the 
range 0.01 < z < 0.16, with a median redshift 0.098. For 
different primary magnitudes, M c , the number of primary 
galaxies and their median redshifts are shown in Table Q] 
For each magnitude bin, the number of primary candidates 
is determined by the interplay between the accessible volume 
given the survey limit and the density of galaxies. The ac- 
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tual number of primaries is further affected by the isolation 
criteria which, for example, tend to reject nearby galaxies 
for which 2-Rinncr subtends a large angle. 

The schematic in Fig. [5] indicates our selection proce- 
dure for potential satellites or "inner galaxies" , and the cor- 
responding selection of the "outer galaxies" used to define 
the background. We assume the satellites of the primary 
galaxy fall within a projected radius, .Rinner (the red circle 
in Fig. [5J). To reduce the background contamination, we ap- 
ply the same cuts in redshift (spectroscopic and photomet- 
ric) as were applied when selecting the primary galaxies, 
but as most of the galaxies within Rinner only have pho- 
tometric redshifts with quite large measurement errors, we 
still cannot distinguish true satellites from projected back- 
ground galaxies. However, the existence of satellites will 
make the number density of galaxies within iti, lncr slightly 
larger than that in the outer blue reference annulus in Fig. [2] 
{Rimer < r < -Router). By counting the difference between 
the number density of galaxies within Rinncr and in the refer- 
ence annulus, we can estimate the number of true satellites. 

An example of the objects we detect around a typical 
primary galaxy is shown in Fig. [3] This image, produced by 
the SDSS finding chart tooQ, illustrates the quality of the 
data and shows that candidate satellites are spatially well 
separated from the light distribution of the primary galaxy. 
The white circle (slightly stretched in this Aitoff projection) 
indicates r = -Rinncr- Within this region we have marked all 
the galaxies in our catalogue with red circles and the subset 
brighter than m r = 20.5, used in our main analysis, with 
yellow boxes. The remaining visible objects within -Rinncr 
are not in our catalogue. Manual inspection with the DR7 
Navigate tool reveals them to be classified as stars. 



3 ESTIMATING THE SATELLITE 
LUMINOSITY FUNCTION 

Once the primary galaxies are defined, their potential satel- 
lites are found from the photometric galaxy catalogue as 
depicted in Fig. [21 For the ith primary galaxy, the num- 
ber of inner galaxies, A^ nner (M), is found by counting all 
neighbouring galaxies within the inner area that satisfy the 
following conditions: at least AMfaint fainter than the pri- 
mary; if they have a spectroscopic redshift, z s , then it should 
satisfy \z c — z s \ < Az s ; or if they only have a photometric 
redshift z p , then it should satisfy \z c — z p \ < a p a*, where 
dp is the error in the photometric redshift as defined below. 
The number of outer galaxies, N° uteT (M) , is determined by 
applying the same conditions to galaxies in the outer area. 
As most satellites of the primary should be projected within 
Rinner of the primary, the number density of inner galaxies 
should typically exceed that of the outer galaxies. The ex- 
cess can be taken as the projected satellite LF of the ith 
primary galaxy, and estimated by 

\ inner 

NT\M) = Nr<*(M) - ±L AC tcr (M), (1) 

where A\ nneY and ^4° utor are the areas of the inner and outer 
regions respectively (excluding sub-regions not within the 

2 http: / / cas.sdss . org/dr7 /en/ tools / chart / chart . asp 
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Figure 4. Estimation of the satellite luminosity function. The 
top panel shows the V-band LF for a single primary galaxy. 
The middle panel shows the mean satellite LF of all primary 
galaxies. The black (thin) and red (dotted) lines give the counts 
of inner and outer galaxies respectively and the blue (thick) 
lines the estimate of the satellite LF. The number of primary 
galaxies contributing to the mean satellite LF in each bin is 
shown in the bottom panel. Here the selection parameters, 
{-Rinner: -Router, AM bin , AM faint , Az s , « p }, are set to the default 
values {-21.25, 0.3 Mpc, 0.6 Mpc, 0.5, 0.5, 0.002, 2.5} 



sky coverage of the SDSS D R7, which we have identified 
using the mask described in iNorberg et all l|201ll )) . 

Because of the survey apparent magnitude limit, we are 
able to probe less of the faint end of the satellite LF for pri- 
maries at higher redshift. To account for this and construct 
an unbiased estimate of the satellite LF averaged over all 
primary galaxies, we count the effective number of primaries 
contributing to each bin of the LF using the weighting func- 
tion 
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Outer galaxies: 

dinner < r p < i?outer and M > M p + AMf aint 

and, if z s available then \z c — z s \ < Az s 
or otherwise \z c — z p \ < a p a* 

Inner galaxies: 

r p < dinner and M > M p + AM faint and, 
if z s available then \z c — z s \ < Az s or 
otherwise \z c — z p \ < a p a* 

Figure 2. Schematic showing the selection of potential satellite galaxies within i?i nn er and of a reference sample within an annulus 
defined by -Ri nnC r < f < Router, used to subtract the residual contaminating background. For both samples we apply the stated redshift 
cuts to reduce background contamination. We also apply the stated absolute magnitude cut to both samples (assuming the neighbouring 
galaxies are at the same redshift as the primary) though this cut is redundant unless Router > 2i?i nncr as otherwise the existence of such 
bright neighbouring galaxies would automatically lead to the exclusion of the primary galaxy. 




Wij(Mj) = 



(M lin 



-Mj) 



Mj < M} im + AMj 

M" m + AM, ' 







Mj > Ml 



lim 



AM, 



where Mj is the central value of each magnitude bin, AMj is 
the half width of the bin, A/" m = m lim -51og 10 (A L )-^(zi), 
Df is the luminosity distance of the ith galaxy and m llm 
is the SDSS galaxy spectroscopic sample magnitude limit. 
For a given primary, the weighting function is unity for all 
magnitude bins in which satellites anywhere in the bin are 
bright enough to be included in the survey. It is zero if all 
satellites within the bin are too faint to be included in the 
survey and ramps between zero and one when only galaxies 
in a fraction of the bin width are accessible to the survey. 
We then define the effective number of primary galaxies, 
ATP™, contributing to the jth bin of the LF as Nf lm = 
^2-Wij(Mj). With this definition, our unbiased estimator 
of the average satellite LF is given by 



N sa -\Mj) 



N piim 



(3) 



In practice, in our study we divide the satellite lumi- 
nosities, Mj, into 20 bins (j — 1,2, •••,20). Furthermore, 
because each primary galaxy in the same bin has a slightly 
different magnitude relative to M c , we choose to show our 
results in terms of the difference in the magnitude of the 
satellite and primary galaxy, AM = M s — M p , which aligns 
the satellite LFs in the same bin. 

The process of estimating the satellite LF for primaries 
in one bin of l^-band absolute magnitude is illustrated in 
Fig. [4] The thin black histogram in the top panel shows the 
number of inner galaxies binned by l/-band magnitude dif- 
ference for one of the primaries. The dotted red histogram 
shows the corresponding number of outer galaxies scaled by 
the ratio of areas A?™ CI /A° utcI . Their difference, which is 
an estimate of the satellite LF in that system, is shown 
by the thick blue histogram. The thin black and dotted 
red histograms in the middle panel show the number of 
inner and (scaled) outer galaxies per primary where the 



number of primaries, Nf 



^2-Wij(Mj), contributing 



at each AM is shown in the lower panel. The heavy blue 
histogram in the middle panel of Fig. [4] shows the estimated 
mean satellite LF for all primaries in the magnitude range 
—21.75 < My < —20.75. The error bars on this mean satel- 
lite LF are estimated by bootstrap resampling of the set of 
primaries. At the faint end of the LF the error bars become 
quite large because of the small number of nearby primaries 
that are able to contribute. If the faintest bin only contains 
one primary then we show a Poisson, rather than the boot- 
strap error. 

For a specific M c , the selection of primaries and counts 
of inner and outer galaxies are determined by the param- 
eter set {dinner, -Router, AMbin, AMf a i n t, Az s , a p }. It is im- 
portant to choose appropriate values for these parameters. 
Here we discuss the physical motivation for our choice of 
parameter values and check that the resulting satellite lu- 
minosity function is robust to reasonable variations in these 
parameters. The various panels in Fig. [5] show the re- 
sults of varying these parameters away from our default 

Choice Of {M c , fi illncr , -Router, AM b in, AMfaint, Az s , Q p } = 

{M c , 0.3 Mpc, 0.6 Mpc, 0.5, 0.5, 0.002, 2.5}. 

The area within which we search for the satellite signal 
is determined by the parameter Rjnncr- For too small a value 
of i?innor, we would lose genuine satellites. Once Ri uaei is suf- 
ficiently large to enclose all the true satellites the resulting 
background-subtracted satellite LF should be independent 
of -Rinner. However, the statistical error in the estimate will 
increase due to increased background contamination. The 
value of 0.3 Mpc is roughly the virial radius of the Milky 
Way, and so this seems a reasonable value to take for the 
dinner of Milky Way-like primary galaxies. One could argue 
for scaling Rj lln cr with the magnitude or type of the primary 
galaxy, but, for simplicity, we set Rj Iln cr = 0.3 Mpc in this 
study except in our parameter tests. In Fig.[5K, we show that 
the effect of varying Rj lln cr between 0.25 and 0.35 Mpc does 
not change the satellite LF significantly. A possible concern 
is that the SDSS data reduction pipeline occasionally mis- 
classifies fragments of the spiral arms of bright galaxies as 
separate galaxies. We have checked that these contaminat- 
ing objects do not make a significant contribution to our 
estimate of the satellite luminosity by excluding all galaxies 
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Figure 3. An example SDSS DR7 image centred on a primary galaxy of magnitude m r = 16.10 at redshift z = 0.074. The white 
ellipse marks r = iJj nnor = 300 kpc. All catalogued galaxies projected within J?i nncr of the primary are marked with red circles. 
Those brighter than our fiducial m r = 20.5 magnitude limit are marked with yellow boxes. The remaining unmarked images 
within -Ri nner are presumed to be classified as stars, which we have verified in this case using the manual SDSS DR7 Navigate 
Tool. 



within 1.5 times the Petrosian R90 radius of the primary 
galaxies. Comparison of the resulting satellite luminosity 
functions shows that they make no significant difference. 

The next parameter, Router, determines the outer refer- 
ence annulus from which we estimate "background" counts. 
An appropriate value for .Router will guarantee a suitably 
local estimate of the backgro und. A local estim ate of the 
background is preferable (see IChen et"afl 120061 ) as galax- 
ies are clustered and, in our case, the mean environment 
of a primary galaxy is also biased by the isolation crite- 
ria that we apply. Fig. [5p shows that, provided the outer 



area is sufficiently large to allow an accurate estimate of the 
background, the resulting satellite LF is robust to changes 
in .Router- We also tested the effect of estimating the back- 
ground using a larger annulus that was disjoint from the 
inner region (from 0.5 Mpc to 0.7 Mpc ) and again found no 
significant difference. 

Besides the physically motivated parameters, we also 
test the parameters of the estimation method. For a spe- 
cific central magnitude, M c , the bin half width, A A/bin, is 
a compromise between having a large enough sample of pri- 
mary galaxies and not distorting the LF due to averaging 
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unlimited rrT m =20.5 



varying surface brightness 
limit (m^ m =20.5) 
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Figure 5. The effect on the estimated satellite LF of varying the parameters {i?i nncr , -Router, AAif bin , AAff aint , Az s , q p } from their default 
values, {—21.25,0.3 Mpc,0.6 Mpc, 0.5, 0.5, 0.002, 2.5}, as indicated in the legends. In addition, panel e shows the effect of changing the 
assumed photometric redshift error from the original cr p to our adopted <r* = max(<r p , 0.05) . Panels h,i and j show the effect of varying 
the apparent magnitude limit of the photometric catalogue and imposing an additional cut on surface brightness (see text for details) . 
In this and in Figs. 6, 7 and 9, the error bars for different datasets have been slightly shifted for clarity. 
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over primaries of differing luminosities. Fig. [5fi shows results 
for a few different AMbin values and indicates that, for our 
choice of binning, the satellite LF by the magnitude differ- 
ence, AAf = M s — M p , any biases are very small. 

The next panel, Fig. [5ji, shows the effect of varying 
the parameter AMfaint, which is important in selecting iso- 
lated primaries. The larger AMfaint, the smaller the number 
of primary galaxies that survive the isolation filter. Hence, 
the value of AMf a i n t is a compromise between avoiding the 
introduction of primary galaxies within groups and gather- 
ing sufficient primary galaxies. We adopt AMfaint = 0.5, but 
Fig- EJi shows that, apart from the truncation of the satellite 
LF brighter than AM = AMf a i nt , the results are, perhaps 
surprisingly, insensitive to changing to AMf a i nt = 0.1 or 1.0. 
To test further the effect of varying the isolation criteria we 
have cross match ed our primary galaxy catalogue with the 
lYang et al.l |2007l ) group catalogue. We find that within the 
DR4 footprint of the Yang et al catalogue only 467 of our 
~ 20 000 primary galaxies for our fiducial value of M c and 
AMf a i n t = 0.5 match with groups of 2 or more galaxies. Ex- 
cluding these group members from our list of primaries has 
essentially no effect on the estimated LF and so we conclude 
that our satellite LF has no significant contamination from 
group members. 

The parameter a p helps us to distinguish genuine satel- 
lite galaxies from background galaxies by excluding galaxies 
that are at a significantly different redshift. If too small a 
value of Qf p is used then we will artificially exclude gen- 
uine satellite galaxies just because the random error in their 
photometric redshift happens to be greater than a p a p . If 
the quoted a p were accurate for all galaxies and the errors 
were Gaussian then a p > 2 ought to be sufficient. However 
the dotted blue and dashed red lines in Fig. [S£ show that 
with both Qf p = 2 and 2.5 the satellite LF is systematically 
underestimated at the faint end. Further investigation has 
revealed that the cause of the sensitivity is that some galax- 
ies with low values of a p in reality have larger redshift errors 
due either to non-Gaussian distributions or inaccuracies in 
(j p . Hence, for our default selection we have been more con- 
servative and set a floor on the photometric redshift error by 
adopting a p — max(<r p , 0.05). Fig. [5f shows that with this 
choice the satellite LF does not depend systematically on 

Qp. 

Fig. [5g shows the dependence of the satellite LF on 
Az s , the maximum allowed spectroscopic redshift difference 
between a satellite and its primary. This should be large 
enough so that satellites are not excluded due to the line- 
of-sight component of their orbital velocities. Our default 
choice is A,z s = 0.002, corresponding to a line-of-sight veloc- 
ity difference of 600 km s _1 . The results are very insensitive 
to this value, mainly as only a small fraction of our potential 
satellites from the photometric catalogue have spectroscopic 
redshifts. 

The final three panels of Fig. [5] illustrate the sensitiv- 
ity of our results to the apparent magnitude and surface 
brightness cuts that we impose on the photometric cata- 
logue. Fig. shows that the satellite LF is systematically 
suppressed at the faint end if all catalogued galaxies are used 
to a faint magnitude limit of my = 21.5, compared to our 
default of 20.5. Brighter cuts also cause some variation but 
in this case the samples are becoming smaller and noisier. 
Figs. [5j and j show the effect of applying cuts in surface 



brightness (mean surface brightness within the Petrosian 
J?50 radius) for two different apparent magnitude limits. For 
a faint magnitude limit of my = 21.5 the faint end of the LF 
is very sensitive to the surface brightness cut. This occurs 
because the catalogue is not complete to my = 21.5 and 
preferentially misses low surface brightness galaxies. With 
the brighter default cut of my = 20.5 this effect is greatly 
reduced (Fig. [5)), indicating much higher completeness and 
little sensitivity to the surface brightness cut. For the r-band 
catalogue, we perform similar tests and find that cuts at sim- 
ilar values to those found for the V-band are appropriate. 

Some of the known Local Group satellites have quite 
low surface brightnesses i|Mated[l998r ) and it is important 
to check that their counterparts would not be missed in our 
analysis by falling below the SDSS detection limit. In Fig. [5] 
we plot the distribution of observed surface brightnesses of 
galaxies around primaries in two different redshift intervals. 
The turnover in these distributions at around Ey « 23 is to 
be expected g iven the intrinsic d istribution of galaxy surface 
brightnesses l|Driver et al-lliooBr i. The distributions for the 
SDSS spectroscopic surve y only become inco mplete around 
Ey ~ 24 mag arcsec" " 2 jstrauss et al.ll200ll ). The surface 
brightness distributions of the subset of Local Group satel- 
lites whose absolute magnitudes are sufficiently bright for 
them to be selected in our catalogue are shown by the blue 
histograms. These can be seen to have surface brightnesses 
that fall near the middle of the measured distribution. 

If the 8 Local Group satellites considered for this study 
were gradually moved to higher redshifts, then only NGC205 
would drop out of our sample by having a surface brightness 
below Ey = 24 mag arcsec -2 before it was lost beneath the 
flux limit. As our sample also includes the SDSS DR7 photo- 
metric subsample, we do actually detect satellites at surface 
brightnesses below that of NGC205, so a conservative esti- 
mate of the incompleteness due to low surface brightness is 
1 in 8. M32 is such a centrally concentrated satellite that it 
would be classified as a star by SDSS, so there is also likely 
to be a comparably small incompleteness at high surface 
brightness in our analysis. 

These combined results show that our method of esti- 
mating the satellite LF is quite robust to changes in the 
parameter values used in the estimation method. There- 
fore We Will USe {dinner, -Router, AM bin , AMf aint , Az s , Q p } = 

{0.3 Mpc, 0.6 Mpc, 0.5, 0.5, 0.002, 2.5} in the rest of the pa- 
per. 



4 RESULTS 

We now explore the dependence of the satellite LF on the 
properties of the primary galaxies. Estimates of the V and r- 
band satellite LF for primaries of magnitude M c = —21, —22 
and —23 are shown in Fig. [7] As the luminosity of the pri- 
mary increases the number of satellites increases at all values 
of AM, and, in addition, the shape of the LF changes. None 
of the luminosity functions are well fit by Schechter func- 
tions, i.e. they are not well described by power laws with 
exponential cutoffs at the bright end. Instead, there is a 
tendency for the LFs to become flatter at the bright end 
and the satellite LFs of the brightest primaries even have a 
local maximum at AAf = 2. Only at the faint end are the 
luminosity functions accurately represented by power laws. 
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Figure 6. A test of surface brightness effects. The black his- 
tograms show the distribution of surface brightnesses, defined as 
the average surface brightness within the half light radius, for 
potential satellite galaxies around primaries at redshift z ss 0.05 
(upper panel) and 2 si 0.1 (lower panel). These are compared 
with the surface brightness distribution (blue histograms) of those 
bright Local Group satellites which would be brighter than our 
my = 20.5 apparent magnitude limit when placed at the redshift 
of the selected primaries. These equiv alent surface b rightnesses, 
computed from the data provided by iMateol Jl998t) , have been 
k-corrected and redshift-dimmed to the redshift of the selected 
primaries. 



We show such fits and list their slopes in Fig. [7] . The vari- 
ety of features in the LFs suggests they will place interesting 
constraints on formation models. 

In Fig. [5] we carry out a comparison of the satellite 
LF of primaries of similar luminosity to the MW and M31 
with data for these two galaxies. The average satellite LF 
of the MW and M31 has often been compared to the oret- 
ical models (e.g. iBenson et al"1l2002l ; ISomervillel |2002| ) and 
used to constrain properties of the model such as the red- 
shift of reionization and the strength of supernova feedback. 
In so doing one implicitly assumes that the satellite LF 
per primary galaxy of the combined MW+M31 system is 
typical of isolated galaxies of similar luminosity. The data 
allow a direct test of this assumption at the bright end, 
My < —14, of the LF. For this comparison, we assume that 
the V-band magnitude s of both the MW and M31 lie in 
the r ange -21.25 ± 0.5 (|Flvnn et al.ll200fj| ; |Gil de Paz et ai] 
120071 ) and compare directly with the average of their V-band 
LFs by plotting on the a;-axis the V-band AM + M c . Over 
the range — 14 > My > — 19 our mean LF has a very sim- 
ilar slope to that of the average of the MW and M31, but 
with almost a factor two fewer satellites at all luminosities. 
Fainter than My = — 14 our estimate becomes noisy due to 



Figure 7. The estimated satellite LFs for different bins of pri- 
mary magnitude as indicated in the legend. The top panel is for 
the V-band and the bottom panel for the r-band. The straight 
lines show power-law fits to the faint ends of the luminosity func- 
tions. Their slopes, a, are given in the legend. Here, the selection 
parameters, {i?mner, Router, A A4"bin, AMf a i n t, Az B , « p }, arc set to 
the default values {0.3 Mpc,0.6 Mpc, 0.5, 0.5, 0.002, 2.5} 



a lack of nearby primaries. The random errors on our es- 
timate of the mean luminosity density are small at bright 
magnitudes, yielding a well-defined estimate of the luminos- 
ity function that provides a very strong constraint on models 
all the way to magnitudes as brigh t as My = —2 . Com par- 
ison with the the oretical models of lBenson et alj (|2002l ) and 
ISomervillel (|2002h highlights the range of predictions. Tuning 
the models to match our new data rather than just the MW 
or M31 may lead to a different assessment of the strength 
of feedback effects in suppressing the formation of satellite 
galaxies. This is particularly apparent when one considers 
the system-to-system variation in the satellite LF. We have 
estimated the intrinsic rms scatter about the mean LF us- 
ing the method detailed in the Appendix. We indicate this 
range with the blue error bars on the cumulative LF in the 
lower panel of Fig.[S]and the mean plus the rms of the differ- 
ential LF by the blue dashed line in the upper panel. Since 
even in the cumulative LF, the mean number of satellites 
per primary is low it is inevitable that the width of the dis- 
tribution includes zero satellites. The wide scatter illustrates 
the danger of just using the MW+M31 to constrain models. 

It is also interesting to see how the satellite luminosity 
function depends on the colour and morphology of the pri- 
mary galaxy. Fig. [9] shows the resulting satellite LFs when 
primaries of l/-band magnitude —21.25 ± 0.5 are split by 
colour and by concentration. In the upper panel we divide 
the primary galaxies into "red" and "blue" subsamples ac- 
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Figure 8. A comparison of the average satellite LF in our sam- 
ple with the satellite LF of the Milky Way and M31. The upper 
panel shows the differential satellite LF of MW-like galaxies. The 
solid line with error bars shows the estimated V-band satellite 
LF of primaries with similar magnitudes to the Milky- Way and 
M31 (M c = -21.25 ±0.5 in the V-band). This is compared to 
the mean LF of the MW and M31 (per central galaxy) in unit 
magnitude bins shown by the r ed points. The best fit power law, 
dN/dM v = 10xl0 01 ( M "+ 5 ). of lKoposov et alJfcOOSh is shown as 
the cyan line. The the oretical predictions of lBenson et al.l (120021 ') 
and ISomervillej (1200211 for 2 rc i on = 10 are shown by the green 
and red lines respectively. The blue dashed line labelled "rms of 
satellites" shows the mean value plus the rms of the LF among 
different primaries. The lower panel shows the same results and 
the observational data in cumulative form. Here, the black error 
bars give the error on the mean cumulative LF while the much 
broader blue error bars indicate the intrinsic rms scatter about 
this mean. 



cording to the well-kno wn colour bimodality i n the colour- 
magn i tude plane (e.g. Strateva et al. l 200 ll; Baldry et alj 
12004 IZehavi et al] 120051 ). Following IZehavi et al l (|2005n T 
we use an equivalent colour cr iterion of °'° (g — r) cut = 
0.19 — 0.24Af r (not identical to IZehavi et al.l as our mag- 
nitudes are K-corrected to z = 0.0 rather than z = 0.1). We 
see that in this bright satellite regime, the LF around blue 
primaries is lower than the LF around red primaries. This 
difference might simply reflect the relative mass of the halos. 
Assuming stellar mass to correlate with halo mass we would 
expect that at a fixed l/-band magnitude blue star forming 
galaxies would be less massive than their red counterparts. 

The lower panel splits the sample into early and late 
type where the early type is defined as having a concentra- 
tion index c ^ 2.6. This division roughly separates early- 
type (E/SO) galaxies fr om late-type (Sa/b/c, Irr) galaxies 
(Shimasaku ct al. 200l]). We see that the satellite LF of late 
types is suppressed with respect to that of the early types. 



Figure 9. The mean satellite LF of different colours (top panel) 
and types (bottom panel) of primary galaxy. The satellite LF of 
early-type or red primary galaxies is shown as a red (dot-dashed) 
line and that of late-type or blue is plotted as a blue (solid) line. 



Given the well known correlation between colour and mor- 
phology this result is consistent with the division by colour. 

We can also use the colour information available in 
SDSS to probe the properties of the satellites. For two bins of 
V-band primary magnitude, Fig. 1 101 shows their satellite lu- 
minosity functions split into red and blue subsamples using 
the same cut in the colour magnitude plane as before. We 
see that at all but the brightest magnitudes the satellites 
are predominately blue and star forming. This is in stark 
contrast with the satellites in groups and clusters where the 
brightest tend to be re d and dead while the faintest are blue 
ijSkibba fc Shethll2009l ). We also note that the LF of the red 
satellites is far from a power law. It has a distinct dip in the 
range from 3.0 < AM < 5.0 and, for the brighter primaries, 
the peak AM « 2.0 that we noted earlier in the total LFs 
is clearly present in the red subsample (and also in the blue 
subsample). 



5 DISCUSSION 

We have constructed a large sample of isolated primary 
galaxies and their fainter neighbours using both the SDSS 
DR7 spectroscopic and photometric galaxy catalogues. The 
samples are sufficiently large that we are able to stack the 
systems and accurately subtract the local background to es- 
timate the mean satellite luminosity function (LF) and its 
dependence on the luminosity, colour and morphology (op- 
tical concentration) of the primary. Our main conclusions 
are: 
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Figure 10. The satellite LF split into contributions from "blue" 
and "red" satellites. The top and bottom panels show the results 
for primaries of V-band magnitude M c = —21.25 and —22.25 
respectively. The "blue" and "red" satellite LFs are plotted as 
blue (solid) and red (dot-dashed) lines respectively. 
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Figure 11. The red line shows th e satellite LF estimated by 
l|Lares. Lambas. &: Do mingucz 2010) for their SO- A-a subsample. 
The solid black and dotted blue lines show our result for a similar 
sample primary galaxies, those brighter than -21.5 in r-band and 
in the redshift range 0.03 < z < 0.1, using AMf a ; nt = 2.0 and 0.5 
respectively. 



(i) The satellite LF is well determined over a range ex- 
tending to approximately 8 magnitudes fainter than the pri- 
mary, for primaries with V magnitudes in the range -20 to 
-23. 

(ii) The satellite LF does not have a Schechter form. After 
a steep decline at the faintest magnitudes, the LF roughly 



follows a fairly flat power law but there is a bump at rela- 
tive magnitude AM ~ 2 which is particularly significant for 
brighter primaries (see Fig. 7). 

(hi) Over the range — 14 > My > —19, the mean satellite 
LF around primaries of My = —21.25 has a similar slope, 
but about a factor of two lower amplitude than the average 
of the combined MW and M31 LFs (see Fig. 8). 

(iv) The amplitude of the satellite LF increases with the 
luminosity of the primary. Over most of the range sampled, 
the increase is approximately a factor of 2 per primary V 
magnitude, but there are significant variations in the shape 
of the function for primaries of different luminosity (see 
Fig. 7). 

(v) The amplitude of the satellite LF also varies with the 
colour and the morphological type of the primary. Red pri- 
maries have more satellites than blue primaries and early- 
type primaries have more satellites that late-type primaries 
(see Fig. 9). 

(vi) Except for the brightest objects, satellite galaxies are 
predominantly blue and star-forming (see Fig. 10). 

As we were completing this work tw o related s tudie s 
were published, both using the SDSS DR7. lLiu et alj (|20ld ) 
used similar selection criteria to ours to construct a sample 
of Milky Way-like primaries and deconvolved for the varia- 
tion of the background to determine the frequency at which 
these Milky- Way like systems host satellites as bright as the 
SMC and LMC. They find that 11.6% host one such satel- 
lite and only 3.5% host two. And they find a mean of 0.29 
satellites per primary. This is in excellent agreement with the 
mean of 0.30 that we find f or satellite s between 2 and 4 mag- 
nitudes (the range used bv lLiu et all ) fainter t han prima ries 
with the magnitude, My = —20.9, adopted bv lLiu et alj For 
the fiducial "Milky Way" luminosity we have adopted here, 
My = -21.25 ± 0.5, we find a slighly larger mean of 0.47 
Magellanic cloud type satellites per primary. 

In a separate study, lLares. Lambas. &: Dominguej 

(20 1 11) estimated cumulative satellite luminosity functions 
and radial density profiles of satellite systems around pri- 
maries brighter than M r — —20.5. When we reproduce the 
selection criteria of one of their samples using our catalogue, 
we find excellent agreement for satellite magnitudes fainter 
than M r = —18.5, but at brighter magnitudes we find a 
significant excess compared to their estimate (see Fig. IA1[) . 
This excess is robust to changes in the value of the isolation 
parameter, AMf a i n t, that we have used. 

The satellite LF probes the smallest scales visible today 
in the hierarchy of galaxy formation. This statistic provides 
a strong test of the ACDM cosmological model, which ro- 
bustly predicts the number of subhalos that could host satel- 
lite galaxies, and a test of galaxy formation theory, which 
determines which of these subhalos are populated by visible 
satellites. The results so far are encouraging. For example, 
the or iginal ACDM galaxy formation model of lBenson et al.l 
(2002) (which predicted the population of ultrafaint satel- 
lites subsequently dis covered in t he SD SS), as well as the 
more recent model of lGuo et alj (I2010T ) predict that bright 
satellites like the LMC and the SMC should be rare. This 
feature appeared to be a shortcomin g of the model whe n 
data were available only for the MW (|Koposov et al.ll2008h . 
The new results for large samples of MW-like galaxies by 
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iLiu et all (|2010h and ourselves now suggest that the MW is 
unusual in having such bright satellites. 

According to standard theory, the satellite LF is es- 
tablished by processes that regulate star formation in small 
halos, namely photoionization of the gas at high redshift 
and supernova feedback, acting on a population of dark 
matter subhalos, itself the result of dynamical evolution 
from a spectrum of primoridal AC DM densi t y pert urba- 
tions. Our analysis and tho s e by I Liu et all (|2010l ) and 
lLares. Lambas. fc Doming^ (|2010l ) reveal features in the 
satellite LF and systematic trends with the properties of the 
central galaxies. These properties encode information about 
galaxy formation processes that will help develop increas- 
ingly refined theoretical models. 
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ies in three equal area regions of the outer annulus are (N bi ), with 
i = 1, 2 or 3. 
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ner area. These two contributions cannot be measured sep- 
arately, but we can estimate the mean number of satellites 

as, 

(N s ) = <A T )-(A b0 ) (A2) 
= <7V T )-/(A out ), (A3) 

where A ou t = Nbi + Ab2 + Ab3 is the total number of back- 
ground galaxies in the outer area and / is the ratio of the 
inner to outer areas. Below we will take / = 1/3 which is 
the case for our default choice of -Router = 2_Ri, mcr . 

Here we are interested in calculating the variance in 
the number of satellites, ((A s — (N B )) 2 ). Starting with equa- 
tion (|A1|) we can write the variance in the total number of 
galaxies as 

<(7V T - (A T )) 2 ) = ((N s + N h0 - ({N s + AM}) 2 } 
= <(iV s - (N B )) 2 ) + ((N h0 - (N h0 )) 2 ) 

+2((N B -(N B ))(N b0 -(N b0 })). (A4) 

If we assume that the number of actual satellites, N B , around 
each primary is uncorrelated with the number of background 
galaxies, Abo, the final cross term vanishes to leave 

<(A S - <iV s )) 2 ) = <(7V T - (AM) 2 ) - {(N b0 - (AW) 2 ). (A5) 

The term, ((A b o — (Abo)) 2 ) cannot be directly measured, 
but, to a good approximation, we would expect it to equal 
the variances, ((N bi — (N bi )) 2 }, of each of the equal area 
portions of the outer annulus. Hence, our final estimate of 
the variance in the number of genuine satellites per primary 
can be written as 

((N B -(N B }f) = ((iV T -(iV T )) 2 )-i^{(iVb I -(A bl )) 2 ).(A6) 

For a given selection of primaries and choice of satellite ab- 
solute magnitude, these terms will depend on the redshift 
of the primary. We find a smooth variation with redshift 
bin and weight the variances according to the contribution 
each redshift makes to the overall estimate of the satellite 
luminosity function to estimate the overall variance on the 
luminosity function. The result is shown in Fig. 8. 



APPENDIX A: ESTIMATE OF THE 
POPULATION VARIANCE 

Here we describe the method by which we estimate the in- 
trinsic variance in the number of satellites per primary. As 
illustrated in Fig. lAll we are only able directly to count the 
total number of galaxies, 

A?t = A s + Abo, (Al) 

where N B is the number of genuine satellites and Abo is 
the number of contaminating background galaxies in the in- 
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